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Aim To compare the reliability of short and long cyto- 
chrome oxidase I gene fragment (COI) in identification of 
forensically important Diptera from Egypt and China. 

Methods We analyzed 50 specimens belonging to 1 8 spe- 
cies. The two investigated markers were amplified by poly- 
merase chain reaction (PCR) followed by direct sequenc- 
ing. Nucleotide sequence divergences were calculated 
using the Kimura two-parameter (K2P) distance model and 
neighbor-joining (NJ) phylogenetic trees. 

Results Although both tested fragments showed an over- 
lap between intra and interspecific variations, long mark- 
er had greater completeness of monophyletic separation 
with high bootstrap support. Moreover, NJ tree based on 
the long fragment clustered species more in accordance 
with their taxonomic classification than that based on the 
short fragment. 

Conclusion In dipterous identification, it is recommended 
to use the long COI marker due to its greater reliability and 
safety. 
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Necrophagous insects can serve as a valuable source of 
information for estimation of minimum post-mortem 
interval (PMI) in legal medicine. Most suitable for foren- 
sic purposes are species from the order Diptera (eg, Cal- 
liphoridae, Muscidae, and Sarcophagidae) (1-4). In PMI 
estimation, an important initial step is correct identifica- 
tion of these insects, which may be difficult by using the 
traditional morphology-based approach (5,6), because 
several forensically important fly species can hardly be 
distinguished morphologically (7-9). The limitations of 
morphological method can be overcome by gene se- 
quences analysis, a fast and accurate method of species 
identification. Molecular analysis requires small tissue 
samples and is relatively insensitive to preservation condi- 
tions (1,10). Different mitochondrial (mt) and nuclear (nu) 
DNA markers are investigated as forensic tools. Howev- 
er, mtDNA is preferred because it can be easily extracted 
even from small or degraded samples (10). In addition, be- 
cause of its strictly maternal inheritance and lack of genet- 
ic recombination, mtDNA haplotype is a good candidate 
for evolutionary and population genetics study. 

Mitochondrial cytochrome c oxidase subunit I (COI) se- 
quences are a rapid and powerful tool for accurate iden- 
tification of species across various taxa (7,1 1 -1 4). Although 
COI has been extensively studied by forensic entomolo- 
gists, resulting in a vast amount of DNA data, there is little 
agreement as to which portion of the gene needs to be 
sequenced. Although the 5'end of COI is also the site of the 
proposed universal animal DNA "barcode" (1 1) and it has 
been successfully used in the identification of many blow- 
fly species (1 2), this approach cannot identify some closely 
related species (12,15). Therefore, to optimize discrimina- 
tion power between closely related species some authors 
suggested multi-gene approach (16,17). Surprisingly, a re- 
cent study using this approach revealed that phylogenetic 
tree based on COI fragment was similar to that based on 3 
different gene fragments (1 6). 

Fragments of the COI sequence that show low sequence 
divergence within species but high divergences among 
species can be employed as taxon "barcodes," and un- 
known samples can be accurately grouped to species with 
reference sequences of the "barcode library" (14,18,19). 
Therefore, it is paramount to evaluate not only discrimina- 
tion power of these COI fragments between closely relat- 
ed species but also between species belonging to more 
than one family, because in a database an unknown sam- 
ple will be compared to all reference samples. In the ab- 
sence of an appropriate reference sample, unknown 



samples will simply group with the most closely matched 
reference sample (20). Thus, it is important to confirm that 
the investigated marker will not only be correctly assigned 
to a species but also that it will be in accordance with the 
traditional morphological classification. Therefore, we eval- 
uated the discrimination power of the short (272-bp) COI 
fragment in the identification of the most forensically rel- 
evant flies (Calliphoridae, Sarcophagidae, and Muscidae) 
originating from Egypt and China in comparison to the 
long (1 1 73-bp) COI fragment, and aimed to gather genetic 
data on common forensically important Diptera. 

MATERIALS AND METHODS 

Samples 

Fifty adult flies belonging to 18 species including 10 spe- 
cies of Calliphoridae, 5 species of Sarcophagidae, and 3 
species of Muscidae were collected during two consecu- 
tive years (1/201 1 to 12/201 2). This study was conducted in 
both Forensic Medicine & Clinical Toxicology department, 
Faculty of Medicine, Suez Canal University, Ismailia, Egypt 
and National Key Laboratory, Basic Medical School, Central 
south University, Changsha, China. All samples were col- 
lected using traps baited with animal remains. Collected 
flies were trapped at different locations in Egypt and China 
(Table 1 ). Samples were identified by entomologists based 
on traditional morphological characteristics (21-25). All 
samples were subsequently stored in 70% ethanol at -20°C 
For comparison, other sequences were retrieved from the 
NCBI database (http://www.ncbi.nlm.nih.gov). 

DNA extraction 

MtDNA was extracted from all samples using Mini Tissue 
Kit (Qiagen, Hilden, Germany) according to the manufac- 
turer's protocol. To avoid possible contamination of fly 
DNA with DNA from ingested proteins and eggs of gut 
parasites, the thoracic muscle of each insect was used as 
the source of DNA, whereas the head and abdomen were 
retained for further analysis. 

PCR amplification and DNA sequencing 

The 272-bp COI gene fragment was amplified using the 
primers 5 '-C AG ATCG AAATTTAAATACTTC-3 ' and 5'-GTAT- 
CAACATCTATTCCTAC-3' and 11 73-bp COI fragment was 
amplified using 5' TACAATTTATCGCCTAAACTTCAGCC 3' 
and 5'CAGCTACTTTATGAGCTTTAGG 3'Details of the prim- 
ers and PCR condition were described in previous studies 
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TABLE 1. Specimens used in the study 



Species 


Code in neighbor-joining tree 


Location 


Accession number 


Chrysomya megacephala (Fabricius, 1794) 


CmC1 


Changsha, China 


KC249623 


KC249673 




CmC2 


Changsha, China 


KC249624 


KC249674 




CmE3 


Ismailia, Egypt 


KC249625 


KC249675 




CmE4 


Ismailia, Egypt 


KC249626 


KC249676 




Cm5 




JX1 87372* 




C. albiceps (Wiedemann, 1819) 


CalbEl 


Alkantra shark, Egypt 


KC249627 


KC249677 




CalbE2 


Alkantra shark, Egypt 


KC249628 


KC249678 




CalbE3 


Alkantra shark, Egypt 


KC249629 


KC249679 




CalbE4 


Ismailia, Egypt 


KC249630 


KC249680 




CalbE5 


Ismailia, Egypt 


KC249631 


KC249681 




Calb6 




AF083657* 




C. rufifacies (Macquart, 1842) 


CrC1 


Changsha, China 


KC249632 


KC249682 




CrC2 


Changsha, China 


KC249633 


KC249683 




Cr3 




JX1 87383* 




C. nigripes (Aubertin, 1932) 


CnCI 


Changsha, China 


KC249634 


KC249684 




CnC2 


Changsha, China 


KC249635 


KC249685 




CnC3 


Guangzhou, China 


KC249636 


KC249686 




CnC4 


Guangzhou, China 


KC749637 


KC74Q687 


Aldrichina graham (Aid rich, 1930) 


AgC1 


Changsha, China 




rv^z^+yooo 




AgC2 


Guangzhou, China 


KC249639 


KC249689 


Lucilia sericata (Meigen, 1826) 


LsCI 


Changsha, China 


[\Lz4yo4U 


vi — >/i OAon 




LsC2 


Changsha, China 


KC 249 641 


KC 2496 91 


L bazini (Seguy, 1934) 


LbCI 


Zhangjiajie China 








LbC2 


Zhangjiajie China 




KT74Q6CR 


L caesar (Linnaeus, 1758) 


LcaCI 


China 


KC249644 


KC249694 




LcaC2 


China 


KC249645 


KC249695 


L. cuprina (Wiedemann, 1830) 


LcuCI 


Changsha China 


KC249646 


KC249696 




LcuC2 


Changsha China 


KC249647 


KC249697 


Lporphyrina (Walker, 1856) 


LpC1 


Changsha China 


KC249648 


KC249698 




LpC2 


Changsha China 


KC249649 


KC249699 


Musca domestica (Linnaeus, 1758) 


MdE1 


Alkantra shark, Egypt 


KC249650 


KC249700 




MdE2 


Ismailia Egypt 


KC249651 


KC249701 




MdE3 


Ismailia Egypt 


KC249652 


KC249702 


M. autumnalis (De Geer, 1776) 


MaC1 


Changsha, China 


KC249653 


KC249703 




MaC2 


Changsha, China 


KC249654 


KC249704 




MaE3 


Ismailia, Egypt 


KC249655 


KC249705 




MaE4 


Alkantra shark, Egypt 


KC249656 


KC249706 




MaE5 


Portsaid, Egypt 


KC249657 


KC249707 


Fanniacanicularis (Linnaeus, 1761) 


FcE1 


Ismailia, Egypt 


KC249658 


KC249708 




FcE2 


Ismailia, Egypt 


KC249659 


KC249709 




FcE3 


Alkantra shark, Egypt 


KC249660 


KC249710 


Sarcophaga albiceps (Meigen, 1826) 


SalbCI 


Changsha China 


KC249661 


KC249711 




SalbC2 


Changsha China 


KC249662 


KC249712 


S. dux (Thompson, 1869) 


SdCI 


Changsha China 


KC249663 


KC249713 




SdC2 


Changsha China 


KC249664 


KC249714 
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TABLE 1. Continued. Specimens used in the study 



Species 


Code in neighbor-joining tree 


Location 


Accession number 


5. Africa (Wiedemann, 1824) 


SaC1 


Changsha, China 


KC249665 


KC249715 




SaC2 


Xining, China 


KC249666 


KC249716 




SaC3 


Changsha, China 


KC249667 


KC249717 




Sa4 




JQ582120* 




5. argyrostoma (Robineau-Desvoidy,1830) 


SargyEI 


Alkantra shark, Egypt 


KC249668 


KC249718 




SargyE2 


Ismailia Egypt 


KC249669 


KC249719 




SargyE3 


Ismailia Egypt 


KC249670 


KC249720 




Sargy4 




JQ582123* 




5. peregrine (Robineau-Desvoidy,1830) 


SperCI 


Changsha, China 


KC249671 


KC249721 




SperC2 


Changsha, China 


KC249672 


KC249722 



*GenBank accession numbers are given for previously published sequences. 



(26,27). Gel electrophoresis was used to isolate PCR prod- 
ucts, which were then purified using QiaQuick PCR Purifi- 
cation Kit (Qiagen, Germantown, MD, USA). Column cycle 
sequencing was performed on both forward and reverse 
strands using ABI PRISM Big Dye Terminator Cycle Se- 
quencing Ready Reaction Kit by ABI PRISM 3730 (Applied 
Biosystems, Foster City, CA, USA) with Big Dye terminator v. 
3.1 as the sequencing agent. 

Sequences analysis and phylogenetic tree construction 

Analysis of DNA sequence variations, nucleotide composi- 
tion, and genetic distances analysis was performed using 
Molecular Evolutionary Genetics Analysis v. 5.10 (MEGA) 
(28). Phylogenetic trees based on the 2 investigated COI se- 
quences were constructed by neighbor-joining (NJ) meth- 
od using Kimura two-parameter (K2P) model implemented 
in the MEGA and tested by 1 000 bootstrap replicates. 

RESULTS 

Both 272-bp and 1 1 73-bp COI fragments were successfully 
sequenced from all 50 insects. The 272-bp and 1 1 73-bp se- 
quences corresponded to positions 2098-2369 and 1513- 
2685, respectively of Drosophila yakuba (GenBank acces- 
sion number X03240). 

Based on 272-bp sequences, 73 were variant and 71 were 
parsimony-informative characters. The nucleotide compo- 
sition showed much higher frequencies of adenine and 
thymine (31 .7% and 37% of total nucleotide compositions, 
respectively) compared with 14.2% of cytosine and 17.1% 
of guanine. NJ analysis was conducted to determine the 
relationships between the analyzed species (Figure 1 and 
Table 2). All species were monophyletic with bootstrap 
support of 99%-100%, except M. autumnalis and F. 



canicularis. Both species could not be separated forming 
one polytypic clade with 61% support. Although Musci- 
dae formed a distinct group with high bootstrap support 
(100%), 272-bp COI marker failed to distinguish between 
Muscidae and Fanniidae. Sarcophagidae family formed a 
distinct group but with low bootstrap support (19%). Cal- 
liphoridae family failed to form a distinct group. At the ge- 
nus level, Lucilia formed a distinct group with 27% support. 
Although Aldrichina grahami belongs to Aldrichina genus, it 
was embedded within Lucilia group. Chrysomya group did 
not join directly with the other group (Lucilia) that belongs 
to the same family. All tested species displayed intraspecific 
variations from 0 to 1 .5% (Table 2). The highest variation was 
observed in C. megacephala and S. africa at 1 .5%. Although 
M. autumnalis samples were collected from 2 countries, 0% 
intraspecific variation was observed. The interspecific varia- 
tions between 18 tested species varied from 1% to 14%. The 
minimum interspecific variations were between M. Domes- 
tica, M. autumnalis, and F. canicularis at 1 %. 

Based on 11 73-bp sequences, 386 were variant and 372 
were parsimony-informative characters. The nucleotide 
composition showed much higher frequencies of adenine 
and thymine (29.9% and 38.8%, of total nucleotide compo- 
sition, respectively), compared with 15.3% of cytosine and 
16.1% of guanine. All tested species were monophyletic 
with full bootstrap supports (Figure 2 and Table 2). Sar- 
cophagidae formed a distinct group with 100% bootstrap 
support. In the Muscidae group, 2 tested families (Musci- 
dae/Fanniidae) could be separated. Calliphoridae family 
failed to form a distinct group. At the genus level, Lucilia 
formed a distinct group with 49% support. Aldrichina gra- 
hami, belonging to the Aldrichina genus, first formed a 
separate group then joined with that of Lucilia with 96% 
support. Chrysomya formed a group with 98% support. In- 
terestingly, Chrysomya groupjoined with Muscidae before 
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joining with other Calliphoridae (Lucilia and Aldrichina). 
All tested species displayed intraspecific variations rang- 
ing from 0 to 2% (Table 2). The highest level was observed 
for 5. africa at 2%. Although samples were collected from 2 
countries, 0% intraspecific variations were observed for M. 
autumnalis. The interspecific variations between 1 8 tested 
species varied from 1% to 15%. The minimum interspecific 
variations were found between L cuprina/L. sericata and M. 
domestica/M. autumnalis/F. canicularis at 1 %. 



lOO. SalhCl 
I&lD)C2 



DISCUSSION 

This study found that although both tested fragments 
showed an overlap between intra and interspecific var- 
iations, long marker showed greater completeness of 
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FIGURE I.The neighbor-joining tree using Kimura's 2-parameter 
model illustrating phylogenetic relationships among 18 fly spe- 
cies based on 272-bp cytochrome oxidase I sequences. Sample 
codes are as in Table 1. Numbers on branches indicate the sup- 
port value. Evolutionary distance divergence scale bar is 0.01 . 



FIGURE 2. The neighbor-joining tree using Kimura's 2-parameter 
model illustrating phylogenetic relationships among 1 8 fly spe- 
cies based on 1 1 73-bp cytochrome oxidase I sequences. Sample 
codes are as in Table 1 . Numbers on branches indicate the sup- 
port value. Evolutionary distance divergence scale bar is 0.01 . 



www.cmj.hr 



24 



FORENSIC SCIENCE 



Croat Med J. 2014;55:19-26 



monophyletic separation with high bootstrap support. 
To our knowledge, this is the first study to provide mo- 
lecular data on forensically important species from Egypt 
and China by using either short 272-bp or long 1 1 73-bp 
fragment of the mt COI gene. The mt COI gene has been 
shown to be a major candidate gene for identification of 
forensically important insects (7,14,27,29). So, before us- 
ing it in real forensic entomology cases, it is worth evalu- 
ating the applicability of different 272-bp and 11 73-bp 
COI genetic markers by using species from the specific 
geographic areas (30). 

As expected, this region of mtDNA had a strong adenine- 
thymine bias, which is characteristic of insect mtDNA 
(6,12). No insertions or deletions were identified within 
the aligned sequences, as was found in studies conduct- 
ed on other mtDNA fragments (6,1 1,31,32). Based on both 
tested COI fragments, C. megacephala and M. autumnalis 
samples were both sequenced from China and Egypt and 
showed minimal variation between populations. However, 
the largest intraspecific variation was observed between 
the species collected from different locations within one 
country. These results are in agreement with the study by 
Harvey et al (20), who tested 1 167-bp COI for identifica- 
tion of Calliphoridae of Australian and South African ori- 
gin. The low intraspecific variation between two countries 



indicates the value of the mtDNA region in interspecific 
distinction (33,34). 

One study suggested that intraspecific variation should be 
<1% and between-species separation >3% (35), where- 
as other studies suggested establishing group-specif- 
ic thresholds (8,1 1). In the present study, results of both 
short and long COI fragments support the idea of estab- 
lishing group-specific thresholds because the 3 investi- 
gated species that belong to Muscidae exhibited the low- 
est interspecific variation, leading to an overlap between 
intraspecific and interspecific nucleotide divergences. In- 
terestingly, although low sequence divergence can result 
in similar haplotypes, which may lead to misidentification 
and a wrong PMI estimate (8), 1 1 73-bp COI was able to dis- 
tinguish between M. autumnalis and F. canicularis without 
bias, but 272-bp COI was not. 

Based on 1 1 73-bp COI gene tree, all species were reciprocally 
monophyletic with full bootstrap support. This observation 
was the same as the analysis based on 272-bp COI fragment, 
except for M. autumnalis and F. canicularis. Surprisingly, trees 
based on both fragments showed that Chrysomya clade did 
not directly join with the other clade belonging to Calliphori- 
dae. This observation may shed light on the importance of 
examining the exact relationship between these groups. 



TABLE 2. Calculated intra- and interspecific divergences expressed as percentage of the analyzed 272-bp (below the di- 
agonal) and 11 73-bp (above the diagonal) cytochrome oxidase I gene fragment using neighbor-joining (NJ) approach with 
Kimura's 2-parameter (K2P) model* 



No 


Species 


N 


VI 


V2 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


1 


C. megacephala 


4 


0-1.5 


0-1.5 




6 


8 


6 


11 


10 


9 


9 


9 


10 


10 


11 


11 


12 


12 


13 


13 


13 


2 


C. albiceps 


5 


0 


0-0.7 


7 




4 


7 


11 


10 


9 


10 


11 


10 


12 


12 


12 


12 


11 


13 


12 


13 


3 


C. rufifacies 


2 


0 


0-1.2 


10 


3 




9 


12 


11 


11 


12 


12 


12 


12 


12 


12 


12 


11 


13 


12 


14 


4 


C. nigripes 


4 


0 


0-0.5 


7 


5 


7 




10 


10 


10 


10 


10 


11 


11 


11 


11 


12 


11 


14 


13 


14 


5 


A.grahami 


2 


0 


0 


12 


11 


11 


11 




9 


8 


8 


8 


9 


14 


14 


14 


12 


12 


14 


13 


13 


6 


L. bazini 


2 


0 


0 


10 


8 


9 


8 


8 




6 


7 


7 


7 


12 


12 


12 


11 


11 


12 


12 


12 


/ 


L. caesar 


2 


0 


0 


8 


8 


8 


10 


9 


7 




6 


6 


5 


11 


11 


10 


11 


12 


11 


12 


12 


8 


L cuprina 


2 


1 


1 


9 


10 


12 


9 


8 


8 


7 




1 


6 


12 


12 


12 


11 


10 


11 


11 


12 


9 


L. sericata 


2 


0 


0 


10 


10 


12 


10 


9 


7 


8 


2 




7 


12 


12 


12 


11 


11 


12 


12 


12 


10 


L porphyrina 


2 


0 


0 


11 


10 


12 


14 


11 


9 


7 


10 


10 




12 


12 


12 


12 


12 


12 


13 


13 


11 


M. autumnalis 


5 


0 


0-0.1 


9 


9 


10 


10 


14 


10 


9 


11 


11 


12 




1 


1 


13 


14 


14 


14 


14 


12 


M. domestica 


3 


0-0.4 


0-0.8 


10 


9 


10 


9 


13 


10 


9 


12 


12 


13 


2 




1 


14 


14 


14 


14 


15 


13 


F. canicularis 


3 


0-0.4 


0-0.6 


10 


8 


9 


9 


13 


10 


9 


11 


12 


13 


1 


1 




13 


13 


13 


14 


14 


14 


S. albiceps 


2 


0 


0 


11 


9 


9 


11 


12 


9 


9 


9 


9 


12 


9 


10 


9 




7 


10 


10 


8 


15 


S. dux 


2 


0 


0 


11 


6 


9 


8 


11 


7 


9 


8 


9 


12 


9 


9 


8 


7 




9 


8 


8 


16 


5. argyrostoma 


3 


0-0.7 


0-1 


12 


12 


12 


12 


11 


9 


9 


11 


12 


12 


10 


10 


9 


9 


6 




10 


9 


17 


S. africa 


3 


0-1.5 


0-2 


12 


10 


9 


9 


10 


10 


9 


11 


12 


12 


12 


12 


12 


9 


7 


9 




9 


18 


S. peregrina 


2 


0 


0 


13 


11 


14 


14 


12 


9 


11 


10 


10 


10 


12 


13 


12 


10 


8 


9 


11 





*Abbreviations: N - number of specimens; V1 - intraspecific variations within 272-bp fragment; V2 - intraspecific variations within 11 73-bp frag- 
ment. 
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Based on 1173-bp COI gene tree, Aldrichina clade pre- 
sented a deviation from traditional taxonomy because this 
species (Calliphorinae) was identified as a sister species to 
Chrysomya ratherthan to Lucilia (16). This pattern of evolu- 
tion was also observed previously based on 28rRNA alone 
(36) and based on COI, CYTB, and ITS2 in a multi-gene ap- 
proach (16). This relation was different from that observed 
based on 272-bp COI, when A grahami was embedded 
within Lucilia tribe. The data obtained by 1 1 73-bp COI phy- 
logenetic analysis were more in accordance with the tradi- 
tional morphological classification than the data obtained 
by 272-bp COI fragment analysis. 

In this preliminary genetic identification of fly species from 
Egypt and China, we found that the long COI fragment 
outperformed the short one in species identification. Since 
the sample size was small, we recommend an evaluation 
of more samples using the same and other loci to confirm 
our findings. In addition, it is important to identify addi- 
tional forensically important fly species and expand such 
analyses to all relevant Egyptian and Chinese species. 
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